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ABSTRACT 

The small grain sizes produced by Type II supernova (SN II) models in young, metal-poor 
galaxies make the appearance of their infrared (IR) spectral energy distribution (SED) quite 
different from that of nearby, older galaxies. To study this effect, we have developed a model 
for the evoluti on of dust content and the IR SED of low-metallicity, extremely young galaxies 
based on Hira shita et alJ (120021) . We find that, even in the intense ultraviolet (UV) radiation 
field of very young galaxies, small silicate grains are subject to stochastic heating resulting in a 
broad temperature distribution and substantial MIR continuum emission. Larger carbonaceous 
grains are in thermal equilibrium at T ~ 50-100 K, and they also contribute to the MIR. We 
present the evolution of SEDs and IR extinction of very young, low-metallicity galaxies. The 
IR extinction curve is also shown. In the first few Myrs, the emission peaks at A ~ 30-50 ^m; 
at later times dust self-absorption decreases the apparent grain temperatures, shifting the bulk 
of the emission into the submillimetre band. We successfully apply the model to the IR SED of 
SBS 0335—052, a low metallicity (1/41 Zq) dwarf galaxy with an unusually strong MIR flux. 
We find the SED, optical properties and extinction of the star forming region to be consistent 
with a very young (age ~ 6.5 x 10 6 yr) and compact (radius ~ 20 pc) starburst. We also 
predict the SED of another extremely low-metallicity galaxy, I Zw 18, for future observational 
tests. We estimate the FIR luminosity of I Zw 18 to be low as £fir ~ 10 Lq, depending 
on the uncertainty of dust mass. Some prospects for future observations are discussed. 

Key words: dust, extinction - galaxies: dwarf - galaxies: ISM - galaxy formation - infrared: 
galaxies 



1 INTRODUCTION 

Galaxy formation is undoubtedly one of the desiderata in modern 
astrophysics, but remains a long-standing unsolved problem. The 
quest for 'primeval galaxies' is still a strong motivation of deep 
surveys by large telescopes and future space missions. By now, 
much effort has gone into the search for galaxies at high redshift 
(z), some of which we believe are experiencing their first episode 
of star formation: Lyma n break objects (e.g., Mada u et alJll996|: 
I Steidel et alll99dfl9 99>: Lyman a emitt ers fe.gi TTLompsoneHrtJ 
1 19951 iHu et alJll998t ISteidel et alJl2000t) : submillimetre galaxies 
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(e.eJHughes et alll998tlEales et alll999ll200(iEiain et all 19991 
Barge r et alJl200CllScott et al]|2002f) : hyper-extremely red objects 
(HER Os) which a re also thought to be strongly enshrouded by dust 
fe.g.. lTotani et all200ll) . 

However, these objects have already produced significant 
amount of stars, and consequently, are chemically enriched by 
heavy elements. Galaxies are thought to be extremely metal-poor 
just after their birth, hence they may have had only a small amount 
of dust. Even in Lyman break galaxies however, there is clear 
evidence that they contain non- negligible amount of dust (e .g., 
lAdelbereer & Steidei2o6ctlSawickl200lllShaDlev et all200ll) . It 
remains difficult to explore the physics of metal-deficient first ob- 
jects in the early Universe because of our observational limitations. 

Among related issues, the extent to which very early star 
formation has been hidden by dust is one of the most interest- 
ing themes. Little doubt remains about the importance of extinc- 
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tion p roperties to evaluate high-z star formation Ce.g.. lSteidel et alj 
1999). Dust grains absorb stellar light and re-emit it in the far in- 
frared (FIR). Indeed, the FIR spectral range represents a unique 
opportunity to study the dust properties and distribution, and sev- 
eral dedicated infrared (IR) space or stratospheric missions are 
planned (ASTRO-F 1 , SIRTF 2 , SOFIA 3 , Herschel Space Observa- 
tory 4 , etc.). 

Furthermore, how dust forms and evolves in primeval galaxies 
needs to be considered in order to understand the chemical and ther- 
modynamical evolution of a metal-poor interstellar medium (ISM). 
Dust grains drastically accelerate the formation rate of molecu- 
lar hydrogen (H2), expecte d to be the most abundant molecule in 
the ISM (e.g., Glover 2003). Hydrogen molecules emit vibrational- 
rotational lines. It is the only coolant in a cool (temperature T <C 
8000 K) dense core of metal-deficient system, since the atomic 
cooling does not work effectively. Hence, the cooling rate of a 
primeval galaxy consequently depends strongly on the dust content 
and H2 abund ance. 

Recently Hiras hita et alj 120021 hereafter H02) have devel- 
oped a model for the evolution of dust content in very young (or 
primeval) galaxies. Dust is expected to be produced at the final 
stages of stellar evolution, and Type II supernovae (SNe II) are 
the dominant source for the product ion of dust grains in young 
star-forming galaxies (e.g.. iDwek & Scalol fl980). Dust formation 
in stellar winds from evolved low-mass stars like red giant branch 
stars (RGBs) or asymptotic giant branch stars (AGBs) can also con- 
tribute considerably, but the cosmic time is not long enough for 
such stars to evolve at high -,z. On the o ther hand, dust is also de- 
stroyed by SN shocks (e.g. Jjones et alll996h : however, H02 have 
shown that the dust destruction is negligible compared to the forma- 
tion rate in young (~ 10 6 -10 8 yr) galaxies. Then they have mod- 
eled the evolution of FIR luminosity and dust temperature in such 
a youn g starburst on the basis of SNe II grain formation model of 
iTodini & Ferrari] <200ll TF01). 

While high-2 primeval galaxies are still beyond the ability of 
present experiments, starbursts in the Local Universe are often as- 
sumed to be good representatives of the star forming activity at 
high redshifts. As such, they are cosmologically interesting as key 
laboratories for studying the ISM, the transport of supernova pro- 
cessed metals, and the chemical evolution of galaxies and of the in- 
tergalactic medium iMadau, Ferrara. & Rees 2001). H02 have paid 
especial attention to the nearby, extremely low-metallicity galaxy 
SBS 0335—052, since this object may be truly experiencing its first 
burst of star formation (age < 10 7 vn lVanzi et al. 2000). The mod- 
els by H02 successfully explain its observed properties such as dust 
mass and temperature, and far infrared (FIR) luminosity. 

In this paper, we further investigate the IR properties of very 
young galaxies, and construct a simple model of mid- to far- 
infrared (MIR and FIR) SEDs of galaxies starting from the H02 
models. As already mentioned above, SNe II are expected to dom- 
inate the dust formation in very young galaxies. The size of dust 
grains formed in SNe II cannot be as large as 0.1-1 fim (TF01), 
and host galaxies are too young for grains to grow in the inter- 
stellar space. TF01 have shown that the size distribution of grains 
produced by SN II is well approximated to be single-sized. A nar- 
row size distribution favoring small grain sizes makes the appear- 
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ance of IR spectral energy distribution (SED) of young galaxies 
drastically different from that of nearby aged galaxies. We, for the 
first time, properly consider the dust size distribution peculiar to 
the very early stage of galaxy evolution, and construct a model of 
the IR SED of very young galaxies. Then, we compare our models 
to the observations in SBS 0335—052, and continue the discussion 
by H02 on evolution of dust content. 

In order to discuss the effects of dust in such galaxies, we must 
first discuss the effect of dust sublimation. Dust grains are expected 
to sublime in H II regions because of the strong radiation from 
the central source (in this case, OB stars). Indeed, the small-sized 
carbon grains (or PAHs) are de ficient in a s trong ultraviole t (UV ) 
radiation environment (e.g.. lBoselli et alll998tlDale et all2001af) . 
However, the radius of the dust cavity caused by the sublimation is 
much smaller than that of the star forming region for classical size 
grains 1 Inoue 2002, and references therein). Further, as will be seen 
later, although our model dust grains are smaller than the classical 
grains, the sublimation process is governed by the energy of each 
injecting photon, because of the small grain cross section. For small 
grains, the peak temperature cannot reach the sublimation temper- 
ature. Hence we can safely neglect the effects of sublimation in our 
discussion. 

The paper is organized as follows: We first model the dust 
content of a young galaxy whose age is less than 10 8 yr in Sj2| 
We show the evolution of the SEDs of young galaxies and their 
IR extinction curve in ^\ Then, we discuss some topics related to 
the infrared properties of young galaxies in f]4] We finally com- 
pare the model predictions with the observed quantities of the ex- 
tremely low-metallicity BCD, SBS 0335-052, in Sgjl] The lowest 
metallicity star-forming galaxy ever known, I Zw 18 (1/50 Zq), 
is another representative of the star forming dwarfs. We present 
a prediction for the SED of I Zw 18 in £14.21 Based on these in- 
vestigations, we estimate how many low-metallicity dwarf galaxies 
and blue compact dwarfs will be detected in future observational 
projects in £14,31 The last section (§|5j is devoted to our conclu- 
sions. 



2 SED MODEL FOR EXTREMELY YOUNG GALAXIES 
2.1 Star formation and dust production 

First we briefly describe the star formation and dust produc- 
tion model proposed by H02. This model treats the evolution of 
dust content in a young galaxy whose age is less than 10 8 yr. 
Hence dust formation in stellar winds from evolved low- mass stars 
like RGBs or AGBs (e.g..lSchmid-Burgk & Scholzll980l:lDeguch: 
198CtlWoodrow & Aumarll982lJGail et alll984lGail & Sedlmav: 



LI ^ 

hi 

vr 



1985U1987i) is negligible because the time interval is too short for 
stars to evolve to such stages. We can therefore safely assume that 
only SN II contributes to the dust formation. For similar reasons, 
the contribution of Type la supernovae to the total supernova rate 
can be also neglected. Further, since dust destruction by SNe may 
be important on a longer timescale of > 10 8 yr 1 Jones et al J 19961 
Hirashita & Ferraral2002l) . we can also safely neglect this process. 
Therefore, the only contributor to the total dust mass (Md us t) in a 
young (< 10 s yr) galaxy is the supply from SNe II. 
The rate of SNe II is given by 



lit) 



) 0(m) dm , 



(1) 
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where ifj(t) is the star formation rate (SFR) at age t (we define 



© 0000 RAS, MNRAS 000, 000-000 



Infrared SED Model of Young Dwarf Galaxies 3 



t = at the beginning of the star formation), <j>(m) is the initial 
mass function (IMF), r m is the lifetime of a star whose mass is 
m, and it is assumed that only stars with m > 8 Mq produce 
SNe II. H02 assume a constant SFR, t[> = tpo, and a Salpeter IMF 
(i/>(m) oc m -2 ' 35 ; the mass range of stars is 0.1-60 Mq) to obtain 
a first estimate 5 . Then the rate of increase of A/d us t is written as 



M, 



dust — ?Tldust7 i 



(2) 



where md ust is the typical dust mass produced in a SN II. We adopt 
m-dust = 0.4 A/q (HF02). We start the calculation from A/dust = 
at t = 0. Because dust production rate is proportional to 7 which is 
proportional to the SFR tpo, Md us t calculated by our model directly 
scales as if>o. The dust mass evolution with ipo = 1 Mq yr -1 is 
shown in Figure 1 of H02. 



2.2 SED construction 

In this subsection, we present the construction of the SEDs for very 
young galaxies. As already discussed above, in their early phase, 
SN II explosions dominate the supply of dust. The galaxies con- 
sidered here are too young to have evolved low-mass stars like 
RGBs or AGBs. These evolved low-mass stars are considered to 
be the source of large-size grains with a power-law like size d is- 
tribution (e.g.. lBiermann & Harwitll98rtfDominik et alll989h . It 
should also be noted that they are too young to make small grains 
grow larger in their diffuse interstellar medium (see, e.g.,,Whittet 
1 19921 pp.223-224). 

We recall that dust grains are schematically divided into sil- 
icate and carbonaceous ones; only these two species are consid- 
ered in this work. TF01 have shown that the sizes of silicate and 
carbonaceous grains formed in SN II ejecta to be about 10 A and 
300 A, respectively. Then the overall size distribution of dust is 
well represented by a superposition of two delta functions, which 
correspond to the sizes of silicate and carbonaceous grains. TF01 
estimated the uncertainty of the grain sizes to be a factor of ~ 2 
within a reasonable range of the relevant thermodynamical param- 
eters. Other sources of uncertainty and their estimates are discussed 
in TF01. Here we should keep in mind the lack of large-size grains 
of ~ 0.1-1 /im in very young galaxies, while a substantial contri- 
bution of such large 'classical' grains is inferred from the observa- 
tion of the Galactic diffuse IR emission or the IR SEDs of nearby 
galaxies. The discrete and small grain sizes make the appearance of 
the SED of young galaxies drastically different from that of aged 
normal galaxies, especially in the IR. We next focus on the most 
important radiation process of small dust grains. 



2.2.1 Stochastic heating of very small grains 

The dust grains originating from SNe II cannot be as large as those 
formed in the atmospher e of evolved AGB stars or grown in the 
diffuse interstellar space iKozasa & Hasegawalll987l ; lKozasa et alj 



5 We assume a relatively smaller mass cut (60 Mq) for the stellar IMF, 
because TF01 only give dust production rate for stars $J 35 Mq. Though 
some assume that the upper mass bound M up = 120 Mq, our frame- 
work is insensitive to the upper stellar cut-off, because the number of such 
high-mass (> 60 Mq) stars is extremely small. The contribution of the 
high-mass stars to UV luminosity is also much smaller than that of the in- 
termediate (3-40 Mq ) mass stars if the stellar luminosity is weighted with 
the Salpeter IMF. In summary, the upper mass limit does not significantly 
affect our conclusions. 



1989; TF01). TF01 clearly showed that the sizes of the silicate 
grains and of the carbonaceous grains are much smaller than that 
of classical large grains. 

Here we note that TF01 used the classical nucleati on the- 
ory fo r the formulation of mass accretion onto grains 1 Fe der et all 
1966). These grain sizes might be even smaller if we use a more 
recent formation theor y of dust g rains from metal-rich gas medium 
(e.g., Ford 1997; Tanaka 2002, and references therein). This re- 
vision may introduce many subtle and perplexing effects in the 
dust formation process, and further investigations are awaited (cf. 
iTanaka et all2002l) . The effect of dust size will be addressed later 
in the context of the origin of near-infrared (NIR) continuum radi- 
ation. 

It is well accepted that very small grains are 'stochastically 
heated', that is they cannot establish thermal equilibrium with 



the ambient radiation fiek 


1 (Purcell 1976 


Aannestad & Kenvon 


1 19791 iDraine & Anderson! 


1985; Guhathakurta & Draine 


1989; 


Siebenmorgen et al. 1992; 


Li&Draine 2001; Draine & Li 2001 



(DL01)): the heat capacity of very small grains is too small to main- 
tain its temperature until the next photon impinges. Hence, the equi- 
librium temperature of a grain cannot be well defined because its 
'temperature' varies violently with time, i.e., it becomes very hot 
each time a UV photon is absorbed by the grain, and subsequently 
cools rapidly to a very low temperature. Consequently, the contin- 
uum shape from very small grains becomes much wider than that 
of a modified blackbody, and the continuum has a strong extension 
to the NIR. 

We adopt this picture and use it to calculate the dust emis- 
sion in the NIR-MIR. Now we formulate the stochastic heating of 
very small dust grains. The heat capacity of a grain has been exten- 
sively studied and modeled by previous papers (DLOh lLi & Drainel 
l200ll and reference therein). We apply the multidimensional De- 
bye model for the heat capacity of a grain, C(T), as introduced by 
DL01. 6 The specific heat of silicate grains is represented by 



C M (T) = (JV atom - 2) ft 



2j - ^ 500 k) + ^ 3 (l500K 



and that of carbonaceous grains is 

T 



C C (T) = (N C - 2)k 



where 



f(_I_\ +2 f( T ) 
^VSeSK/ n \ 2504 k; 



/»(*) 

fUx) 



V" dy 



dx 



Q exp(y/x) 



(3) 
(4) 

(5) 
(6) 



and k is the Boltzmann's constant (the subscripts 'siT and 'C rep- 
resent the silicate and carbonaceous grains, respectively). Namely, 
as for silicate grains, two-thirds of the vibrational modes are dis- 
tributed according to a Debye model with n = 2 and Debye tem- 
perature = 500 K, and one-third of the modes are described by 
a Debye model with n = 3 and S = 1500 K. For carbonaceous 
grains Nc — 2 out-of-plane modes and 2(Ac — 2) in-plane modes 
are set. This model excellently reproduces C(T) of each grain 
species (see Figure 1 of DL01). We adopt p c = 2.26 [g cm~ 3 l 
toraine & Leell984 and p aU = 3.5 [gem -3 ] <Li & DrainehoOlh . 
i.e., Nc = 1.14 x 10 23 [cm" 3 ] and A at0 m = 8.5 x 10 22 [cm" 3 ], 



6 As the authors have already pointed out in their subsequent paper, there 
is a typographic error in their Equation (10) in DL01. 
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Wavelength X [/^m] Wavelength X [fim] Wavelength X [fim] 




1 10 100 1 10 100 1 10 100 

Wavelength A [fim] Wavelength X [/im] Wavelength X [yitm] 



Figure 1. Evolution of IR SED of a very young galaxy. The size of the star forming region is rgp = 30 pc. Solid and dotted lines represent the SED history 
for a galaxy with the size of silicate grain a s n = 10 A and 6 A, respectively. For Figures 0and|2] constant star formation rate of SFR = 1 Mq is assumed. 



respectively (AT a tom is the number density of atoms). Here we iden- 
tify bulk graphite grain with our carbonaceous grain. 

The (effective) equilibrium temperature, T C q is computed by 
equating grain heating by photon absorption and cooling by photon 
emission (see Equations (6.1) and (6.2) o f lDraine&Leell984l) : 



/ hc\ 4 

u\cQ(a, X)dX = 60u I — x 
XtykJ 



Q(a,\)\- 



exp [he/ (AfcT cq )] 



■dX. 



(7) 



where u\ is the interstellar radiation field energy density, a is a 
grain size, and h, a and c are Planck's constant, Stefan-Boltzmann 
constant, and the speed of light, respectively. We adopted the values 
presented in lDraine & Le for the absorption efficiencies, 

Qahs(X). In practice, the above equation can be approximated by 
setting Q at UV (on the lhs of Equation^ to be unity and Q at FIR 
(on the rhs) to behave like Qfir oc A~ 2 . Then we obtain 



and 



M5u 



TvkJ \ 960iY(2-irAa)hc 



1/6 



u\dX , Qfir 



2nAa 



(8) 



(9) 



Here we adopt j4 s ii = 1.34 x 10 3 [cm] for silicate grains 



iDrapatz & Michellll977t) and A c = 3. 20 x 10~ 3 [cml for car- 
bonaceous grains (scaled to fit Q m (X) o flDraine & Leel 19841) . 

The actual dust temperature T has a very wide distribution 
around T cq . As we will see later, while carbonaceous grains can 
be in equilibrium with the ambient radiation field, silicate grains 
cannot establish the equilibrium, even in the intense radiation field 
(high photon density) of a young galaxy (J3|and §3). 



2.2.2 Emission 



Following Draine & Andersort ll985t) . with this C(T), we calcu- 
lated the temperature distribution of dust as a function of size a 
with Monte Carlo simulations for incident UV photons. 

We calculated u\ from the history of OB star luminosity of 
H02 and the size of the star forming region vsf- We use the OB 
star luminosity calculated in H02 as the total luminosity of stars 
with m > 3 Mq. We start from rsF = 30 pc as the size of a star- 
forming region in young galaxies (H02), but other sizes are also 
examined later. 

We approximate d the sp ectrum of the UV field by the formula 
of lDraine & Bertoldl ll996h for simplicity and scale it with lumi- 
nosity by factor x(t) as 



«a = 



25 
" 6 



10 J x X 
1 cm 



25 
2 



10 5 x A 
1 cm 



13 

3 



xlO 9 x(i) [erg cm 3 cm 



(10) 
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Figure 2. Same as Figure^but the size of the star forming region is rgp = 100 pc. 



The UV radiation field strength of t he solar neighborhood, m nbh , 
is u NBii = 4.0 x lCT 14 [erg cm -3 ] iDraine & Bertoldll99fih . The 
intensity factor \(t) is determined as a function of the galaxy age 
from the OB luminosity history and tsf (Eauation ll 1> 7 . 

For the geometry, we assume that the OB stars are concen- 
trated in the centre of the system, and dust is assumed to be dis- 
tributed as a shell, so as to express the fact that dust absorbs the UV 
photons very effectively in star forming regions iHunt et all200ll) . 
Moreover, as seen later, such a "simplest" geometry is successful 
enough to explain the SED of SBS 0335—052. We regard the radius 
of the sphere, r as the size of the spherical star forming region, tsf- 
Then, the flux is calculated by J(t) = LoB(i)/(47rr| F ). Then the 
UV energy density incident on dust grains, u, is 



= X{t)u 

c 



(11) 



The spectral UV energy density u\ is calculated from the assumed 
spectrum of the UV radiation field (Eauation ll Ot . 

Such a shell-like geometry of dust distribution ma y result from 
the drift induced by the radiation pressure. Gail & Sedlmavr 1 1979) 
show that grains in H II regions can be accelerated nearly to ~ 



7 We note that, for a high- z galaxy, we should include CMB in the radiation 
field. Small grains spend most of their time with their lowest temperature, 
which is determined by CMB. This slightly affect the submillimetre SED 
of high-z galaxies. An extensive discussion on the heating of dust by CMB 
is given in Ferrara et al. 1 1999a). 



10 km s . The crossing time i C ross can be estimated as 



3 x 10° 



30 pc 



10 km s" 



[yr]- 



(12) 



Therefore, the dust grains can be swept into a shell surrounding the 
star-forming region within a timescale o f 10 6 yr. Stellar winds can 
also sweep dust grains (see Inoue 2002, and references therein for 
the creation of central dust cavity). We will go back to this issue in 
line with the application of our model to the local low-metallicity 
star forming dwarf galaxies. 

The rate at which a grain absorbs a photon with energy E ~ 
E + dE is expressed as 



d 2 P r, i ,x 2 A 
dEdt= Q * bs{a ' X)na UX Wc 



(13) 



he 

T 



47T 

T 



The heating is represented as follows: 

l-T 

C{T')dT' , (14) 

'T 

where T is the peak temperature achieved by a grain hit by a photon 
with energy hc/X, and To is the grain temperature just before ab- 
sorption. If more photons are absorbed by a grain simultaneously, 
the lhs will become the total energy of the photons. On the other 
hand, the grain cools through radiation as 



dT _ 3tt 
dt ~ ~ aC(T) J Q 



3 (o,A) B x {T)d\. 



(15) 



From Equations d 1 31 . 1 141 . and <15> . we obtain the sample path of 
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a grain temperature as a function of time. A sufficiently large set 
of Monte Carlo re alizations yields the dust temperature distribu- 
tion l Draine & Anderson 1985). Though there could be faster algo- 
rithms I Guhathakurta & Draine 1989), we adopt this direct method 
in order to examine the cooling behavior and timescale of very 
small dust grains in an extremely intense radiation field of a young 
galaxy. 

The total mass of each grain component is given by TF01. 
The mass ratio we adopt here is M s n : Mrj = 0.56 : 0.44 (cf. 
H02). With this value and material density of each species (p s n = 
3.50 gem -3 and pc = 2.26 gem -3 ), we obtain total numbers of 
grains, Ni 



Ni 



3Mdu S t fi 



(16) 



where subscript i denotes the species of dust, sil or C, and /j is the 
mass fraction of dust of species i. 

Monochromatic luminosity of a grain with radius a and 
species i, Lf^ ln '' (X), is, then, obtained by 



dPi(x,<h) 



dT 



dT, (17) 



where dPi(x, a i)/dT is the temperature probability density func- 
tion. In the present case, at is fixed for each dust component to a 
specific value in the following range: a s a — 6-10 A = 6-10 x 
10~ 8 cm and a c = 200-300 A = 2-3 x 10~ 6 cm, i.e., the dust 
size distribution function is well approximated by a Dirac's delta 
function 8 for each species, i, 

drii 
da 



(a) = UiS(a — at) ■ 
Total IR luminosity from species i, L\ R , is 



(18) 



L} R = iV l / igT^AJdA 



AirafNin / / Ql bs (A)B A (T) 



dPj(x, Oj 

dT 



dXdT, (19) 



and we have imposed the constraint Lm — £\ L\ R ^ Los{t); 
if Lm exceed Lob we set Lm, = Lob - Final continuum IR emis- 
sion from a young galaxy is calculated by superposing each of the 
continuum from silicate and carbonaceous grains as a function of 
starburst age, t, 

Lm,x{t)=J2 L ™W 

i 

= ^4™?iV i7r J Ql hs WBx(T) dPt( * T ai) dT , (20) 

i 

where in the rhs starburst age is incorporated with Ni and 
dPi(x, a l )/dT through M dust {t) and x(t). 



2.2.3 Extinction 

If the dust opacity is very large, then self-absorption occurs, and 
even MIR radiation from dust is absorbed by the dust itself. In 



8 The results will not be affected by the detailed choice of the functional 
form: e.g., Gaussian with small standard deviation, instead of delta func- 
tion, changes the results little. The discrete nature of the size distribution is 
fundamental here, i.e., slowly varying functions such as a power-law yield 
significantly different results in contrast. 



this work we did not try to solve the detailed radiative transfer but 
rather treat the extinction as given by a one-zone screen mo del. The 
resulting extinction curve is different from that of Draine &Led 
1 1984) because of the different dust size distribution. In this case 
dust opacity and extinction are expressed as 



Tdust(A) 7raiQ abs ,i(A) mrsF ■ 



I(X) = I (X)e 



-TdustM 



(21) 



(22) 



The absorbed light is re-emitted at longer wavelengths, mainly in 
the submillimetre, and consequently, the SED is deformed by the 
self-absorption. The final SED is obtained via this self-absorption- 
re-emission process with an obvious constraint of the energy emit- 
ted in the IR, Lir, which must be conserved. 

The screen geometry might seem an oversimplification, but 
many observations support screen-like dust configuration in in- 
tense starbursts (e.g.. lMeurer et al]|l995l) . Moreover, as discussed 
in EI2.2.2I dust can be swept into a shell surrounding the star- 
forming region. The possible dumpiness of dust may also increase 
the dus t opacity t\ of the galaxy (for a discussion see lFerrara et alJ 
1999b), as well as it may allow the UV photons to leak out of the 
galaxy. Although such an anisotropy is difficult to model, it is re- 
markable that the simple geometry we adopted can reproduce a 
number of the observed properties of low metallicity star forming 
dwarfs as we will see later. This aspect will be further considered 
in j4"T2l 



3 RESULTS 

3.1 Infrared SED evolution 

We present the results of the evolution of IR SEDs for a very young 
galaxy with the age of 10 6 -10 8 yr in FiguresQand|2| To illustrate 
the results, and facilitate comparison with H02, we concentrate on 
the value SFR = 1 Mq yr -1 . We set the silicate dust grain size to 
Osil = 10 A and 6 A (solid and dotted lines, respectively), because 
grain size can be changed by a factor of about two depending on 
the conditions of grain formation. The size of carbonaceous grains 
can also be altered for the same reason, but it affects the SED very 
little because the grains are in a thermal equilibrium. Thus, we only 
show the results for ac = 300 A in FiguresQand|2| In the present 
model, the radiation field is typically 10 4-5 times stronger than in 
the solar neighborhood. 

In the case of smaller silicate dust size, continuum radiation 
produced by stochastically heated grain extends toward NIR wave- 
lengths. This effect is quite sensitive to the grain size of a < 10 A, 
in our case that of silicate grains. In contrast, the so-called 'very 
small grains (VSGs)' in the local normal (viz. solar metal abun- 
dance) star forming galaxies are generally thought to consist of 
carbonaceous grains or their fragments. In the early, optically thin 
phase of a young galaxy, these small silicate grains produce a 
prominent band feature around 9.7 fim in emission. It is worth 
mentioning that even in the extremely intense radiation field of 
X = a few x 10 4 -10 5 , such a small grain cannot establish ther- 
mal equilibrium with the ambient radiation, because the cooling 
time of a grain (typically t coo i ~ 10~ 2 s at T ~ 500 K) is much 
shorter than the mean interval of photon injection onto the grain 
(At ~ 10-10 3 s). 

IPlante & Sauvagel 120021) suggested that the absence of 
unidentified infrared band (UIB) features in the intense UV envi- 
ronment of SBS 0335—052 implies that stochastically heated small 
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Figure 3. Left panel: Evolution of IR extinction curve of a very young galaxy. Size of the star forming region rgp = 30 pc. Right panel: Same as the left 
panel but rgp = 100 pc (note that the ordinate is an order of magnitude different from that in the left panel). 



grains are unimportant in such galaxies. However, as we have ex- 
plained, small silicate grains considered here may be uncorrelated 
with UIBs. Small silicate grain is naturally expected in SNe (TF01), 
so we should properly take into account the effect of stochastic 
heating of small grains in starburst systems. Indeed, until the ex- 
tinction becomes significant, their contribution is dominant in N- 
MIR wavelengths in our results. 

We also stress that the dust temperature of carbonaceous 
grains is very high, T > 100 K in the early phase. Because of the 
lack of classical size (~ 0.1-1 /im) large grains, submillimetre con- 
tinuum radiation is weak, and the peaks are located around the MIR 
during the first a few tens of Myrs. Such a high dust temperature is 
also p redicted for giant forming ellipticals by iTotani & Takeuchl 



curve does not change in time. Hence, strictly speaking, Figure [5] 
represents the evolution of dust opacity, rather than the evolution of 
extinction curve. 

By our definition, A\ is independent of the grain size di, since, 
from Equations 1 16i and 12 1 1 . A\ is obtained as 



A x = 1.086r dust (A) oc ( M^Il) 9^1 (23 ) 



3.2 Opacity evolution 

After the early phase characterized by the radiation from hot dust 
(age < 10 7 yr), opacity of dust becomes large, and consequently 
self absorption of dust gets quickly severe in the N-MIR regime. 
Finally, the observed dust temperature decreases to 50 K or less. 
In particular, the 9.7 fim silicate band causes a prominent absorp- 
tion feature even for a young age t ~ 10 7 yr. The absorbed light 
is then re-emitted in the FIR-submm wavelengths, which results in 
a cooler effective temperature. Consequently, the temperature his- 
tory is consistent with that of H02 which was derived simply from 
energy considerations. 

We show the evolution of the IR extinction curves for the cases 
rsF = 30 pc and 100 pc in Figure|3| We should note that, since dust 
is supplied only from SNe II (i.e., the dust size distribution does not 
vary) within the timescale considered here, shape of the extinction 



and Qabs.i scales as Q a b s ,i oc a^. Hence, results are the same for 
di = 6 A and 10 A. 

Here we should keep i n mind that the larg e opacity makes 
radiation pressure effective jFerraraetall ll991). and may cause 
the expansion of the star-forming complex and eventually, the on- 
set of galactic wind at some stage in the galaxy evolution. Then 
the opacity will become smaller than shown in Figure [5] In order 
to take it into account, we should construct a consistent frame- 
work which can treat dynamical and chemical evolution as well 
as radiation. We regard dynamical modeling as the next stage and 
do not go further in this work. Some comprehensive discussions 
on galactic winds in the similar condition can be found in, e.g., 
IPlante & Sauvaeel l2002t) . More general properti es of the wind in 
formin g dwarf galaxies are discussed in, e.g., Ferr ara & Tolstovl 
(2000). 
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Figure 4. The model SED of an extremely young galaxy calculated by our model and the observed SED of SBS 0335—052. Red curve represents the model 
SED. Filled squares are the observed values, while open squares are the upper limits. Open triang les in the near infrared are the net dust contributions to 
the total flux (see main text). Blue curve is the ISOCAM spectrum obtained bv lThuaneta?]<1999l) . Left panel: All photometric values are converted to the 
monochromatic luminosity vL v [Lq]. Right panel: Same as the left panel but expressed by observed flux densities Si, [/iJy]. 



4 DISCUSSION 

4.1 A starbursting dwarf SBS 0335-052 

4.1.1 Obseri'edSED 

Here we adopt our model to the local BCD, SBS 0335-052 with 
its unusual IR SED and strong flux at MIR in spite of its low metal- 
licity (1 /41 Z Q ). Since lHunt eTailfeOOll) estimated the SFR of this 
galaxy to be about 1.7Moyr _1 from Bra line intensity, we scale 
the OB star luminosity of H02 with this value. The result is shown 
in Figure |4| The photometric data presented in filled and open 
squares in Figure [31 are mainly taken from Table 1 of iDale et alJ 
1 2001b). The points at 10.8 and 21 /im obtained by Gemini/OSCIR 
and the data (and more strict upper limits) at 60, 6 5, and 100 umob- 
tained by ISOPHOT are from Table 1 of lPlante & Sauyasd <2002l) 
The IS OCAM spectrum (blue curve) is taken from Thuan et al. 
1 1999). In the NIR, we should correct the contribution of stellar 
light and nebular continuum. By assuming the stellar light fractions 
to be ~ 37 % (2 /im) a nd ~ 6 % (4 nm) an d applying the nebular 
emission coefficients of lov & Lester 1 1988), Hunt et al. 1 2001) ob- 
tained the net dust contribution of ~ 13 % and ~ 67 % at 2 /im and 
4 /tm, respectively. We show these values by open triangles in Fig- 
ure|4| In the left panel, the data are converted to the monochromatic 
luminosity vL v [Lq], while in the right panel we show in the unit 
of the flux density in Figure|4] 

From the NIR spectrum, th e age of th e stellar population is 
estimated to be < 5 Myr IVanzi et al 12 000). In order to reproduce 
an extreme NIR-FIR SED of SBS 0335-052 within the strict age 
limit, stellar energy radiated in the optical wavelengths must be 
effectively converted to FIR radiation. It suggests a very large dust 
opacity, requiring a very compact SF region (rsF = 20 pc). 

The small size of the SF region is also consistent with the ob- 
served high dust temperature ~ 80 K. SBS 0335—052 has a strong 



N-M IR continuum without UIB or PAH features iThuan et alJ 
I1999J) . In our framework, this is consistent with a small silicate 
grain size a sl i = 6 A. Further, the MIR peak and lack of a 
prominent submillimetre component strongly support the absence 
of classical-size large grains. The size of the carbonaceous grains 
is suggested t o be ac = 200 A in order t o reproduce the ISO PHOT 
observation of lPlante & Sauvaed <2002h . |Plante & Sauvagel <2002h 
also reached to a similar conclusion for the lack of classical size 
grains from a different radiative transfer model. 

A significant absorption feature at A ~ 10 /im suggests a 
large dust opacity. It is a natural consequence of the large IR lu- 
minosity of this galaxy, and is also consistent with visual extinc- 
tion of Ay = 12-20 mag for the embedded star forming region 
fThuan et all 19991 : iHunt et all200ll) . This aspect will be discussed 
in fl4~L2l 

As we mentioned in £12.2.21 the configuration of dust may be 
shell-like within a timescale of ~ 10 6 yr. The estimated age of the 
starburst of SBS 0335—052 is comparable to the crossing timescale 
of the system. 

With these physical values and under our assumptions, we 
have succeeded in reproducing the overall SED of SBS 0335—052 
with the starburst age t ~ 6.5 Myr. Thus, all the peculiar proper- 
ties of SBS 0335—052 are consistent with those of a genuine young 
galaxy in the Local Universe. This success strongly supports the va- 
lidity of our model, despite its simplicity. 



4.1.2 Extinction curve 

IHunt et all 1200 ll) have shown that the extinction properties of 
SBS 0335—052 can be well described by the extinction curve of 
iLutJ ll999l) . advocated for the extinction of the Galactic Centre 
(GC), assuming a simplifying assumption of screen-like dust ge- 
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Figure 5. The IR extincti on curve obt ained from our model, with age t ~ 
6.5 Myr (s olid line). ThelLutJ 1 1999ft extinction curve, scaled to Ay = 
16.3 mag fauntet alJl200ll) is also shown (dashed line). Open symbols 
represent the extinction deduced semi-empirically from the 'bare' model 
SED and the observed data: squares are the semi-empirical extinction curve 
derived from the photometry, and triangles are from the ISOCAM spectrum. 

ometry. In their evaluation of extinction, Hunt et al. assumed a 
superposition of two modified blackbodies for the continuum ra- 
diatio n, only the c ooler of which was exti nguished by the screen. 
iThuan et alj 1 19991) and lDale et alj 1200 lbl) also assumed modified 
blackbody for the unextincted continuum to explore the properties 
of the IR radiation of this galaxy. 

However, as we have already seen in the previous section, we 
cannot expect equilibrium between the UV radiation field and very 
small dust grains, which are thought to be the main source of the 
continuum radiation in the N-MIR. This is also pointed out by 
iDale et alj |2001b). Bearing this in mind, here we reconsider the 
IR extinction properties. 

The IR extinction curve predicted by our model is presented in 
Figure[5] This curve is the one of a galaxy with the same parameters 
used in Figure|4] 

First, it is roughly consistent with the lLutJ 119991) extinction 
curve , scaled by a visual extinction Ay — 16.3 mag iHunt et ail 
l200ll) . Indeed, the GC curve is quantitatively well reproduced by 
our model with SFR = 1.7M© yr _1 , the star forming region size 
rsF ~ 20 pc, and the above starburst age. Since our model pro- 
vides the extinction curve only at wavelengths A > 1 fim, if we 
extrapolate the model extinctio n of if -band to F-band by the ratio 
tabulated in Rieke & Lebofskv 1 1985), we obtain the visual extinc- 
tion strength Ay ~ 18 mag. It is als o equal to the value calcu- 
lated by the relation A K = O.lLAv ICardelli et ai]ll989l) . This 
value is in good agreement with the estimate of Hu nt et al., al- 
though lower than the value Av = 30 mag by Plan te & Sauvaed 
1 2002). This may be caused by the difference of the radiation mod- 
els used for calculation. Here we should keep in mind that the Ay 
measured at the optically visible region of SBS 0335—052 is small 
(Ay ~ 0.55 mag; e.g.. llzotov et alJll997l) . typical of dwarf irreg- 



ulars. So, the peculiar characteristics of SBS 0335—052 may be 
governed by the heavily embedded star-forming region. 

Though the model extinction curve is roughly consistent with 
the GC one in the NIR regime, our model predicts a flatter slope 
than the GC one in the extrapolation to the o ptical. Our ext i nction 
curve in the NIR regime has a slope similar to lCardelli et alj ll989l) 
(oc A~ 16 ). lHuntetai]l200ll) found that the extrapolation of Lutz's 
GC curve to t he optical by the slope of oc A -1,75 , close to that of 
ISeatorl 119791) (see the extension of the dashed line to the optical in 
Figure|5J. The flatter slope of our extinction curve is a natural con- 
sequence of the fact that the continuum radiation is not produced 
from the dust in equilibrium with ambient UV field. The shape of 
the non-equilibrium continuum is much broader than that of modi- 
fied blackbody radiation. The intensity of our continuum radiation 
in the N-MIR is stronger than that of the modified blackbody if 
we fix the intensity at about 20 /im, hence, we need a stronger ex- 
tinction (flatter slope) at N-MIR wavelength range. The stronger 
extinction at the 9.7 /im broad silicate feature is also explained in a 
similar way. Such behaviour is required to reproduce the absorption 
feature in the observed SED of SBS 0335-052. 

In order to examine these issues in more detail, we estimated 
the extinction with our 'bare' model SED, calculated without tak- 
ing into account the self-absorption, rather than using the modi- 
fied blackbody continuum. The deduced extinction is represented 
by open squares and triangles in Figure [5] We call the newly de- 
rived extinction the 'semi-empirical extinction curve'. Unfortu- 
nately, since the NIR fluxes from dust emission contain a large 
uncertainty depending on the nebular and stellar contribution, we 
cannot strongly constrain the shape of the curve, but the NIR part 
of our model extinction curve shows rough agreement with the 
semi-empirical estimates. Around the 9.7 /im silicate feature, the 
model extinction curve shows very good agreement with the semi- 
empirical one. This result vindicates the performance of our model, 
as well as the simplifying assumption of a screen-like dust geome- 
try. 

However, around the 18 /im silicate feature, our model sig- 
nificantly underestimates the extinction. Correspondingly, our 
model SED overproduces the flux around 20 /im, compared 
with the observation of Plante & Sauvage (Figure |4j. Many 
IR galaxie s show the strong a bsorption feature at 18 /im (e.g., 
Dudley & Wvnn-Williams 1997). The discrepancy stems from the 
shape of the extinction curve around 20 /im. We cannot reproduce 
such a fe ature within our prese nt framework based on the dust emis- 
sivity of Draine & Lee 1 1984). This riddle will be investigated in 
future work. 

Except the above discrepancy, our model extinction curv e nat- 
urally explains most important properties of the GC curve of lLutzl 
1 1999). The success of our model in reproducing the extinction 
curve strongly supports the composition of dust supplied by SNe 
II as predicted by TF01. 

4.2 Prediction for I Zw 18 

I Zw 18 has long been a record holder as the lowest metallicity star- 
forming galaxy in the Local Universe (Z = 1/50 Zq), as well as a 
textbook example of a dwarf starburst. Hence, in the context of our 
model, it is an intriguing object to disc uss as another representative 
of extremely low-metallicity galaxies. Ostlin (2000) estimated its 
distance to be 12.6 Mpc, and we use this value for our discussion. 
Whether I Zw 18 is really a newly formed system or not has been 
a matter of debate. The discovery of AGB stars may 
suggest an underlying old population, but its contribution is almost 
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Figure 6. The predicted SED for I Zw 18. The NIR observed data are taken from Hunt et al. 1 2003). Upper limits are calculated from the limits of/RAS. Left 
panel: The same as Figure|4] the photometric values are converted to the monochromatic luminosity vL v [Lq]. Right panel: Same as Figure|6] but expressed 
by observed flux densities S„ [p.3y] . 



certainly not dominant <Hunt et alliool) . For this reason, we there- 
fore regard I Zw 18 as a forming galaxy for the first approximation. 

Recently , fro m HST WFPC2 narrowband imaging, 
ICannonetai1j2002h derived various important physical properties 
of this galaxy. They showed that the 08 ionizing equivalent number 
of this galaxy is 2886 for the total galaxy. It corresponds to a total 
OB star luminosity, Lob — 4.9 x 10 s Lq. From the Ha image, 
we estimate the total size of the SF regions of rsF ~ 50-100 pc. 
The SFR of I Zw 18 is ~ 0.04 M Ft yr " 1 , and for some patches, 
Av ~ 0.5 mag iCannon et al]l2002h . We use these values to 
constrain the model prediction. 

From the measured exti nction, I Zw 18 ma y also contain 
a significant amount of dust. ICannon et al] l2002h estimated the 
total dust mass Afd us t — 2-5 x 10 3 Mq. This value is con- 
sistent with the dust-to-gas ratio expe cted from the metallicity 
iLisenfeld & Ferrardll998l : ICannon et all2002h . The dust mass of 
our model reaches this value at the starburst age t ~ 10 7-7 ' 5 yr, an 
age consistent with the observationally suggested major SF age in 
I Zw 18 (e.g.. lMartirll996r.lHunt et all2003[) . Here, we note that a 
smaller size of the SF region makes Ay exceed 0.5 mag and is not 
permitted. 

We show the resulting SED in Figure |6] The photometric 
data in the NIR are taken from lHunt et alj |2003f). For the case of 
I Zw 18, interstellar conditions are quite different from those in 
SBS 0335—052, in spite of their similar metallicities. Since the vi- 
sual extinction Av is very low even in its star forming patches, the 
radiation reprocessing from UV to FIR is inefficient. It naturally 
leads to a low IR luminosity Lir, < 10 7 ' 5 Lq. For the same rea- 
son, the silicate band features are observed in emission in the SED 
of IZw 18. 

It is worthwhile to focus on the differences between I Zw 18 
and SBS 0335—052. Indeed, these two representative low- 
metallicity star-forming dwarfs are significantly different in spite 



of their almost same metallicities. First, the SFR of I Zw 18 is 
much smaller than that of SBS 0335—052. This means that the 
dust production rate of I Zw 18 is much smaller than that of 
SBS 0335—052. In spite of their optical similarity, I Zw 18 has 
rather normal SFR while SBS 0335—052 has extremely vigorous 
star formation. Such an active/passive star formation variety may 
play an important role in understanding their activity I Hu nt et alJ 
2002). Second, the apparent SF region size of I Zw 18 is much 
larger than that of SBS 0335—052. This difference could be pro- 
duced by the dust sweeping by the radiation pressure from the cen- 
tral stellar cluster. When the age of the system becomes older, the 
dust layer may be swept far away from the star-forming complex, 
and the opacity may get smaller. Since I Zw 18 has an older age 
than SBS 0335—052, the dust layer of the I Zw 18 can be trans- 
ported further away. The small dust opacity of I Zw 18 might be 
understood by the combination of the above two effects. 

Inspection of Figure |6| tells us that the peak of the SED of 
I Zw 18 is located around 100 /im both in the most optimistic 
and pessimistic predictions. However, the detection in these wave- 
lengths is limited by Galactic cirrus. A shorter wavelength (~ 
60 /im) observation will be dramatically easier. The MIR observa- 
tion from space seems exciting, and crucial to examine the validity 
of our prediction. In the NIR, dust contribution may be negligible 
in the NIR. The fraction of dust emission appears to be smaller 
than that of SBS 0335—052. This suggests that either of stellar and 
nebular contribution is dominant in the NIR. 

We discuss the possibility and prospects for future observation 
in the next subsection. 



4.3 Prospects for future observations 

As we have seen in the above, even in the intense UV radiation field 
of very young galaxies, small silicate grains suffer from stochas- 
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tic heating and consequently have a broad temperature distribu- 
tion. Hence, silicate grains mainly contribute to the MIR continuum 
emission. Carbonaceous grains have larger sizes and can be in ther- 
mal equilibrium, but the strong UV field makes their equilibrium 
temperature very high (~ 50-100 K), and they also contribute to 
the MIR. These facts indicate that both (restframe) MIR and FIR 
surveys are necessary to examine the dust amount of galaxies over 
wide range of metallicity. 



Our models show that the observed flux of I Zw 18 is expected 
to be 10-15 mjy at 100 u,m. Thi s agree s very well with a simple 
argument bv lKamava & HirashitalkoOlh . who showed that the ex- 
pected flux for I Zw 18 is 7-100 mjy at 100 Lim. Typical detection 
limits for the forthcoming FIR surveys are 1-100 mjy at ~ 80 Lira 
and 10-100 mjy at ~ 150 /im. We expect ~ several x 10 2 detec- 
tions of such local low metallicity galaxies by some space observa- 
tions in the near future such as the SIRTF SWIRE Legacy program 
and ASTRO-F FIR all-sky survey, if we assume the local luminos- 
ity functions for star-forming dwarf galaxies. However we should 
be aware of the cirrus confusion limit, which makes it difficult to 
detect IR extragalactic sources at wavelengths A > 100 /im. In or- 
der to detect these galaxies at A = 100-200 /im, a large aperture 
(<f> ~ 2-4 m) FIR space telescope like Herschel is required to ob- 
tain a small diffraction limit to overcome the cirrus confusion lim- 
its. Herschel Space Observatory will be one of the optimal instru- 
ments for this purpose (on the estimation of the FIR confusion limit 
for various forthcoming instruments, see e.g.. llshii et alJl2002t) . In 
any case, such a sample will be useful for the systematic studies 
on 'primeval galaxies', through the investi gation of other po ssibly 
'young' local dwarfs (e.g., Hopk ins et alj|2002t ICorbin & Vacc3 



Furthermore, our result for hot dust also offers a potentially 
important cosmological insight into high-z galaxies: if dust in high- 
z object is as cool as usually assumed, then we cannot expect a high 
comoving SFR density because in a realistic cosmology the cosmic 
IR background (CIRB) spectrum strongly constrains the high-z IR 
emission (Takeu chi et alJ l2001a). However, very hot dust in such 
high-z object reconciles the CIRB constraint with high SFR, and 
allows vigorous star formation hidden by dust in the early Uni- 
verse iTotani & Takeuchi 120 021. Some studies attempt to predict 
the number of submillimetre galaxies at high-z (z = 2-3) from 
both empirical and theoretical viewpoints (e.g.. | G uiderdoni et al 
I l998t iTan et alJI 19991: iMaeliocchetti et alj|200lt iTakeuchi etaf 
2001b). Most models assume galaxy SEDs similar to those of local 
dusty giant starbursts like M82 or Arp 220. However, we advocate 
our theoretical framework for very young galaxies to make it possi- 
ble to treat dwarf starbursts, because present galaxy formation mod- 
els suggest a vast number of dwarf galaxies in the early Universe. 
Hiras hita & Ferrari] ^2002) provided a comprehensive framework 
of very high-z (or, say, primeval: z > 5) dwarfs in the submillime- 
tre and millimetre, and showed the ALMA detectability. They also 
suggest the existence of dust warmer than the typical Galactic en- 
vironment because of the intense radiation field. This emphasizes 
the importance of our SED modeling which is applicable to young 
galaxies with intense interstellar radiation fields. With our SED 
modeling, thus, the advent of ALMA will enable us to compare 
our models with observations of such low mass, bona ride primeval 
galaxies. 



5 CONCLUSIONS 

We have constructed a simple model of mid- to far-infrared (MIR 
and FIR) SEDs of galaxies to investigate the IR properties of 
very young galaxies , based on the galaxy evolution model of 
Hirashita et all feOOl H02). 

Small grains cannot establish thermal equilibrium with the 
ambient radiation field because of their very small heat capacity. 
We find that, even in the intense UV radiation field of very young 
galaxies, small silicate grains still cannot be in equilibrium, and 
consequently have a broad temperature distribution. Hence, sili- 
cate grains contribute mainly to the MIR continuum emission. Car- 
bonaceous grains are larger and can be in thermal equilibrium, but 
the strong UV field makes their equilibrium temperature very high 
(~ 50-100 K); hence they also contribute to the MIR. 

We calculated the SEDs of very young, low-metallicity galax- 
ies as a function of starburst age. In the first several Myrs, such a 
galaxy shows very hot dust emission with a peak at A ~ 30-50 /im, 
because the self absorption is not very pronounced in this phase. 
Later, the SED changes its appearance to a more common shape, 
dominated by cooler dust emission, caused by the strong self- 
absorption of MIR light and re-emission at submillimetre wave- 
lengths. The evolution of the dust opacity and the shape of the IR 
extinction curve were also given. 

We succeeded in reproducing the MIR-FIR SED of 
SBS 0335-052, a low-metallicity (1/41 Z Q ) galaxy with an un- 
usual IR SED and a strong MIR flux. SBS 0335-052 is consistent 
with a very young starburst of age 6.5 x 10 6 yr and compact size of 
~ 20 pc. The model results of young age and compact size are also 
consistent with the optical properties of the galaxy. We also suc- 
cessfully reproduced the extinction curve of SBS 0335—052, which 
agrees well with the semi-empirical extinction curve deduced from 
the observational data and our unextincted 'bare' continuum SED. 
We obtain the visual extinction strength Av ~ 18 mag, consis- 
tently with Hunt et al.'s estimate of Av = 16.3 mag. 

For future observations, we then made a prediction for the 
SED of another extremely low-metallicity galaxy, I Zw 18. As in 
the case of the optically visible part of SBS 0335—052, I Zw 18 
has a low visual extinction, which suggests a small dust opacity. 
Hence it follows that the FIR luminosity of I Zw 18 is probably 
low as LpiR ~ 10 7-7 ' 5 Lq, depending on the observational con- 
straint of dust mass. We expect the peak of the IR SED of I Zw 18 to 
be around 100 ^im, but Galactic cirrus may hamper its detection at 
these wavelengths. Shorter wavelength observations will be much 
more effective. 

We conclude that both of MIR and FIR surveys are necessary 
to infer the dust amount of galaxies with a wide range of metallic- 
ities. About several hundreds of local low metallicity galaxies may 
be detected in the such as the SIRTF SWIRE Legacy program and 
ASTRO-FFIR all-sky survey, and even more detailed observations 
will be possible by Herschel and other 2-3 m class FIR space tele- 
scopes. Such samples will be useful for the systematic studies on 
the local populations of very young galaxies, as well as for their 
counterparts at high-z. 

In addition, redshifted peaks of very hot dwarf galaxies will 
be found in the submillimetre wavelengths (~ 50(1 + z) /im, and 
z > 5), and we conclude that the advent of ALMA will enable us 
to observe such low mass, bona Gde primeval galaxies. 
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